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Abstract 

As gravitational lensing measurements become increasingly precise, it becomes nec- 
essary to include ever higher order effects in the theoretical calculations. Here we 
show how the difference between the shear and the reduced shear manifest them- 
selves in a number of commonly used measures of shear power. If we are to reap 
the science rewards of future, high precision measurements of cosmic shear we will 
need to include this effect in our theoretical predictions. 

Key words: Cosmology, Lensing, Large-Scale structures 
PACS: 98.65.Dx, 98.80.Es, 98.70.Vc 



1 Introduction 

Weak gravitational lensing by large-scale structure (Mellier, 1999; Bartelmann 
& Schneider, 2001) is becoming a central means of constraining our cosmolog- 
ical model (see e.g. Hoekstra, Yee, & Gladders, 2002; van Waerbeke, Mellier 
& Hoekstra, 2005, for the current status). Continuing effort has yielded in- 
creasingly stringent control of systematic errors and ever larger surveys are 
decreasing the statistical errors rapidly. As the precision with which the mea- 
surements are made increases, one demands ever higher fidelity in the data- 
analysis, modeling and theoretical interpretation. 

One such area is the computation of the theoretical predictions for various 
statistics involving the measurable 'reduced shear' (Mellier, 1999; Bartelmann 
& Schneider, 2001) 
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g = - for \g\ < 1 (1) 



E-mail: mwhite@berkeley.edu 



Preprint submitted to Elsevier Science 



2 February 2008 



where 7 is the 'plain' shear and k is the convergence which is related to the 
projected mass along the line of sight. Both k and 7 are defined through the 
Jacobian of the mapping between the source and image planes. If the lensing 
is weak then \k\ <C 1, and to lowest order (7 = 7 and the power spectra of the 
shear and convergence are equal: Cj = C*. Thus measurements are usually 
compared to predictions of Cf derived from models of the non-linear clustering 
of matter. This relation has corrections however, which need to be taken into 
account if precise comparisons with data are to be made in the future. While 
there has been previous analytic work on this subject (Schneider et al., 1998; 
Dodelson & Zhang, 2005) and some numerical work on magnification statistics 
(Takada & Hamana, 2003; Barber & Taylor, 2003; Menard et al., 2003), there 
has been no comprehensive investigation of this effect in simulations. In this 
short paper we present some results, derived from simulations similar to those 
described in (White & Vale, 2004), on the size of the corrections for a variety 
of commonly used statistics on the angular scales where most current lensing 
work has focused. 



2 Simulation 



The results are derived from 7 N-body simulations of a ACDM cosmology run 
with the TreePM code (White, 2002). Each simulation employed 384 3 equal 
mass dark matter particles in a periodic cubical box of side 200 /i~ 1 Mpc run 
to z = with phase space data dumped every 50 h^Mpc starting at z — 4. 
The cosmology was the same for each simulation (f2 mat = 0.28, Vl B h 2 = 0.024, 
h = 0.7, n — 1 and as = 0.9) but different random number seeds were used 
to generate the initial conditions. The purpose of these runs was to allow a 
study of the statistical distributions of lensing observables. Each box was used 
to generate 16 approximately independent lensing maps, each 3° x 3°, using 
the multi-plane ray-tracing code described in (Vale & White, 2003). Each map 
was produced with 2048 2 pixels and then downsampled to 1024 2 pixels 2 . Thus 
a total of 112 maps or 1008 square degrees were simulated, all with the same 
cosmology. We report on the results for sources at z ~ 1 (specifically all at 
a comoving distance of 2400 h~ 1 Mpc) here. Other data can be obtained from 
http: / / mwhite.berkeley.edu/Lensing. 



2 Tests with maps downsampled to 512 2 pixels indicate that our results are con- 
verged at the several percent level in the ratios that we quote, with better agreement 
at larger scales as expected. 
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Fig. 1. (Left) The smoothed shear variance, Var[|g|^], as a function of smoothing 
scale R. (Right) The relative difference of the smoothed shear variance computed 
for plain shear and reduced shear, as a function of smoothing scale R. The points 
show the mean difference, averaged over 112 maps, each 3° x 3°, while the error 
bars indicate the standard deviation. 

3 Results 



3. 1 Two point statistics 



The lowest order information on large-scale structure comes from studying the 
2-point statistics of the shear field. A number of different 2-point statistics have 
been used in the literature, each with their own strengths and weaknesses. We 
survey the most commonly used statistics here. 

We first show the results on the variance of the shear, smoothed on a range 
of angular scales. Though it has numerous drawbacks in interpretation, this 
was one of the first statistics computed from shear maps due to its ease of 
computation. The left panel of Figure 1 shows the (reduced) shear variance 
as a function of smoothing scale. For each point we compute a map of the 
amplitude of the reduced shear, smooth the map with a 2D boxcar of side 
length R and then compute the variance of the resulting map. The points 
are highly correlated. Since we have full information about the shear and 
convergence from the simulation we are able to compare this with the plain 
shear results which are usually computed. The right panel of Figure 1 shows 
the difference of this statistic computed for the reduced shear (g) to that 
for plain shear (7). We see that the reduced shear variance is larger than 
the usually predicted shear variance by a non-trivial amount. The bias very 
gradually drops as we go to larger scales (not shown here) but we are unable 
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to follow it to convergence due to the finite size of the fields simulated. 

Naively one might imagine that the corrections would be smaller than shown 
in Figure 1 because k is small 'on large scales'. However k is only small when 
averaged over a sizeable region of sky, and this smoothing does not commute 
with the division in Eq. (1). This fact also means that perturbative calculations 
must be used with care, because the convergence or shear amplitude can be 
quite large on small scales leading to a break down of the approximation. 

Though the smoothed shear variance has the worst behavior of the two point 
functions we consider, we will see qualitatively similar behavior below. In 
general this tendency for the small-scale power to be increased over the the- 
oretical prediction for plain shear needs to be considered before attributing 
"excess small-scale power" to additional physical effects e.g. intrinsic galaxy 
alignments. A signature of the effect might be the difference in the sensitivity 
of the various statistics that we compute here. 

Another popular measure of the power is the aperture mass (Schneider et al., 
1998). This is a scalar quantity which can be derived from an integral over 
the shear 

M ap (# ; R) = J d 2 e Q{\e\-R) lT {e + e ) (2) 

where 7t is the tangential shear as measured from 9 and Q is a kernel. We 
have chosen the £ = 1 form for definiteness, so 

Q(r = X /R;R) = r 2 (l - r 2 ) for r < 1 (3) 

and Q vanishes for r > 1 . The results for the aperture mass variance are shown 
in Figure 2. Note that the correction is quite large on small scales, but drops 
rapidly on scales above a few arcminutes. 

Finally we show results for the 2-point correlation function of the shear in 
Figure 3. This is the best behaved 2-point statistic that we consider, since it 
explicitly eliminates the contribution from small angular scales 3 . We define 
£ = (7+7+) where 7+ is (minus) the component 71 of the shear in the rotated 
frame whose x-axis is the separation vector between the two points being 
correlated. The component at 45° is called 7 X . We find that the correction to 
the correlation function is quite small on the scales shown here. The correction 
to the other correlation function, (7x7x), is even smaller than for £. 

It is not too surprising that £ receives smaller corrections than Var[M ap ] when 
we recall that M ap probes smaller scales than R by a factor of roughly 3 
(Schneider et al., 1998). A smaller piece of the difference is that M ap can be 

3 For this reason the 2-point correlation function is less sensitive to small-scale 
observational systematics than other quantities such as Var[M ap ]. 
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Fig. 2. The relative difference of the M ap variance computed for plain shear and 
reduced shear. The points show the mean difference, averaged over 112 maps, each 
3° x 3°, while the error bars indicate the standard deviation of the ratio. 



expressed as an integral over £ extending all the way to zero lag 4 , making M ap 
slightly more sensitive to the difference between g and 7 on small scales. A 
similar argument holds for Var[|<7|#], with the effect being more pronounced. 
A reduction in the sensitivity to small scale power and to shot-noise could be 
obtained from differencing these measures between scales, with an associated 
loss of power. 

These are all of the 2-point functions commonly derived from shear data - 
the power spectrum has only been derived by two groups (Pen et al., 2003; 
Brown et al., 2003) despite being almost ubiquitous in forecasts of the po- 
tential of future lensing experiments (e.g. Hui, 1999; Huterer, 2001; Benabed 
& Bernardeau, 2001; Hu, 2001; Weinberg & Kamionkowski, 2002; Munshi & 



Wang, 2003; Refregier et al., 2004; Takada & White, 2004; Huterer & Takada, 



2005). 



4 The variance of any convolution of k must have weight at zero lag. 
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Fig. 3. The relative difference of £(r) computed for plain shear and reduced shear. 
The points show the mean difference, averaged over 112 maps, each 3° x 3°, while 
the error bars indicate the standard deviation of the difference. The effect on (7x7x) 
is smaller and is omitted. 

3.2 Three point statistics 

Since the lensing maps are non-Gaussian there is information beyond 5 the 
2-point statistics. We show in Figure 4 the dominant 3-point function for 
equilateral triangles as a function of scale. The 3-point function is (7+7+7+) 
with 7+ denned with respect to a line joining the triangle center to each 
vertex (Takada & Jain, 2002; Zaldarriaga & Scoccimarro, 2003; Schneider & 
Lombardi, 2003). In the weakly non-Gaussian limit one can argue that the 
corrections to the 3-point function from considering reduced shear could be 
very large. Counting powers of our small parameter, k or 7, the lowest order 
contribution vanishes in the Gaussian limit. Thus the first no n- vanishing term 
is of order k a (Schneider et al., 1998; Dodelson & Zhang, 2005). The first 
correction from the reduced shear also comes in at order /t 4 . However we see 
that the difference is not significantly larger than in the case of the 2-point 
functions on the scales shown here. 



5 This information is not, however, independent. We find a strong correlation be- 
tween the amplitude of the 2- and 3-point functions on scales of several arcminutes 
in our maps. 
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Fig. 4. The relative difference of the 3-point correlation function, C+++> f° r equilat- 
eral triangle configurations computed for plain shear and reduced shear. The points 
show the mean difference, averaged over 32 maps, each 3° x 3°, while the error bars 
indicate the standard deviation. 

We also show, in Figure 5, how the configuration dependence of the 3-point 
function is modified. The 3-point function for triangles with sides Q\ = 2' 
and 9 2 = 3' is plotted as a function of the cosine of the included angle. We 
have chosen here a triangle with no special symmetries in order to illustrate a 
"typical" case. Note that the correction is dependent on angle, showing that 
the shear and reduced shear have different shape dependence which should be 
taken into account in cosmological inferences (e.g. Ho & White, 2004; Dolney, 
Jain & Takada, 2004). 



4 Conclusions 

Deflection of light rays by gravitational potentials along the line of sight in- 
troduces a mapping between the source and image plane. The Jacobian of 
this mapping defines the shear and convergence as a function of position on 
the sky. In the absence of size or magnification information neither the shear 
nor the convergence is observable, rather the combination g — 7/(1 — k) is. 
Since on small scales k can be non- negligible this introduces complications in 
predicting the observables of weak lensing. We have illustrated specifically the 
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Fig. 5. The relative difference of the 3-point correlation function for one triangle 
size, £+++(#1 = 2', 62 = 3'), computed for plain shear and reduced shear. The points 
show the mean difference, averaged over 32 maps, each 3° x 3°, while the error bars 
indicate the standard deviation. 

effect of using g rather than 7 on a number of commonly used statistics of 
weak lensing. We found that the correlation function is least affected, as ex- 
pected since it specifically eliminates small-scale information. The variance of 
the shear amplitude, smoothed on a scale R, is the most dramatically affected, 
and the amplitude of the effect declines very slowly with increasing smoothing 
scale. 



Existing lensing experiments derive most of their cosmological constraints from 
angular scales of a few arcminutes to a few tens of arcminutes. The corrections 
described here are below the current statistical errors, and so should not affect 
existing constraints. The next generation of experiments may have sufficient 
control of systematic errors, and sufficient statistics, that this effect will need 
to be included in the analysis. Note that even for the correlation function the 
difference between the commonly predicted plain shear result and the reduced 
shear value is comparable to the level of systematic control that we need to 
achieve to study dark energy! If we are to reap the science rewards of future, 
high precision measurements of cosmic shear we will need to be able to include 
this effect in our theoretical predictions. 

M.W. thanks M. Takada for comparison of results on the shear 2-point func- 
tion and D. Huterer, B. Jain and P. Zhang for helpful comments. The simula- 
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tions used here were performed on the IBM-SP at NERSC. This research was 
supported by the NSF and NASA. 



References 

Barber A.J., Taylor A.N., 2003, MNRAS, 344, 789 
Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 291 
Benabed K., Bernardeau F., 2001, Phys. Rev., D64, 083501 
Brown M.L., et al., 2003, MNRAS, 341, 100 
Dodelson S., Zhang P., 2005, preprint [astro-ph/0501063] 
Dolney D., Jain B., Takada M., 2004, MNRAS, 352, 1019 
Ho S., White M., 2004, Astrophys. J., 607, 40 

Hoekstra H., Yee H.K.C., Gladders M.D., 2002, New Astron. Rev., 46, 767 

Hu W., 2001, 2001, Phys. Rev., D66, 3515 

Hui L., 1999, Astrophys. J., 519, 9 

Huterer D., 2001, Phys. Rev., D65, 063001 

Huterer D., Takada M., 2005, preprint [astro-ph/0412142] 

Mellier Y., 1999, Ann. Rev. Astron. Astrophys., 37, 127 

Menard B., Hamana T., Bartelmann M., Yoshida N., 2003, A&A, 403, 817 

Munshi D., Wang Y., 2003, Astrophys. J., 583, 566 

Pen U.-L., et al., 2003, MNRAS, 346, 994 

Refregier A. et al., 2004, Astron. J., 127, 3102 

Schneider P., Lombardi M., 2003, A&A, 397, 809 

Schneider P., van Waerbeke L., Jain B., Kruse G., 1998, MNRAS, 296, 873 

Takada M., Hamana T., 2003, MNRAS, 346, 949 

Takada M., Jain B., MNRAS, 344, 857 

Takada M., White M., 2004, Astrophys. J., 601, LI 

Vale C, White M., 2003, Astrophys. J., 592, 699 

van Waerbeke L., Mellier Y., Hoekstra H., 2005, A&A, 429, 75 

Weinberg N., Kamionkowski M., 2002, preprint [astro-ph/0210134] 

White M., 2002, Astrophys. J. Supp., 143, 241 

White M., Vale C, 2004, Astroparticle Physics, 22, 19 

Zaldarriaga M., Scoccimarro R., 2003, Astrophys. J., 584, 559 



9 



